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SYNOPSIS 


A Thesis fiatitled "lacompressible Slip Plows with 
Permeable Boundaries" submitted m Partial Jhlfilment 
of the ’lequiremeats for the Decree of Doctor of 
Philosophy by S V Sx-GHIDMiUDA to the Departmeit of 
Mathematics, Indian Institute of Te/^hnology Lanpur, 

August 1975 under the Supervision of Dr R K Jam 

Incompressible flow of fluids over porous boundaries have 
attained a nev/ dimension with the worics of Beavers and Joseph 
Before their work on flows over poirous boundaries, research workers 
tacitly assumed that there was no relative motion of the fluid with 
reference to the boundary The purpose of having a porous boundary 
was either to have a controlled mjection of fluid or suction of 
fluid Beavers and Joseph (l967) proposed an adhoc boundary 
condition at a permeable wall a linear relation bet veer the velocity 
and Its gradient at a permeable boundarj;^ Theoretical piedictions 
on flow rates based on this model were expernner tally confirmed 
Theoretical justifications of the model were advanced by Taylor (l97l) 
and Saffman (l97l) aJid subsequent experimental work by other research 
workers have confirmed these observations 

Por flows m ducts and channels with permeable walls, there 
are some natural queries which comes to one's mmd Por instance, 
will the slip condition in the two regions - entry and fully 
developed region - be the same*? How entry length will be affected 
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and what will be the flow in the entiy region The present 
thesis addresses itself to some of these problems 

Aq adhoc variant to the condition of Beavers and Joseph is 
proposed indicating a variable slip from tne leading edge of a 
peimeable flat plate As a consequence slip at the wall is a 
point dCunction in entiy region of the channel Tnis variant has 
an alternative theoretical basis A search for similar solutions 
over a permeable flat plate necessitates slip to vary directly as 
the square root of the distance from the leading edge 

This thesis is divided into six chapters Chapter I gives 
motivation and description of the proposed study A brief review 
of wort done in related topics relevant to the discussion in the 
thesis is given in Chapter II 

Incompressible flov/ of a viscous fluid over a permeable flat 
plate IS studied in Chapter HI As in the Blasius problem, a 
search for similar solutions is made This requires tnat the slip 
velociiy at the porous wall be proportional to the square root of 
the distance from the leading edge along the plate, the proportiona- 
lity constant depending on the permeability of the plate The 
method of solution adopted is the numerical integration The 
velocity profiles are determined for different values of a certain 
non-^imensional number which depends on the properties of the porous 
material These profiles lie between the Blasius profile and 
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the ouxve f' = 1 , where is the aon-^imeasiorial velocity 
It IS observed that the velocity p2X)files are sharper than the 
Blasius profile Mother method of solution of the same 
problem is also suggested Under simple transformation of the 
dependent and independent variables, the governing differential 
equation is reduced to another boundaiy value problem which is 
easier to handle In this foimulation one need not take lacoursc 
to a method of searching the initial conditions as for example is 
done for the Blasius problem itself 

Chapter IV contains a discussion of the entry flow of an 
incompressible viscous fluid in a porous walled channel Ihe 
governing equation of motion are the Prandtl boundaiy layer 
equations, and here ag a in allowance for slip at the wall is made 
The method of solution i§ tne method of von Karmau-Iohlhausen The 
resulting differential equation has been solved for the two cases 
when the velocity profiles m the entry region of tne cbaiinel are 
si mi lar and are nan-similar. The entry length corresponding to 
two definitions are determined - one through the boundary layer 
and the other through the fully developed velocity profile In 
the former case, the entry length is that length vhen the boundary 
layer attains half the width of the channel In the latter it is 
that length when the velocity in the channel attains the value of 
the velocity of the fully developed flow at the centre of the channel 
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For both the cases, the entjy length is shovra. to reduce with a 
non -dimensional parameter which is proportional to the permeability 

Chapter 7 has a discussion of the Benard problem of a differ- 
entially heated peimeable block Ihe problem here is to determine 
the linear instability limit as a function of the properties of the 
porous material The normal mode analysis has been made and the 
resulting system of equations have been solwed using the Galerkin 
method The critical Eayleigh numbers so determined are considerably 
reduced for the porous case as compared with napeimeable walled case 

The final chapter gives some comments on the thesis content 



CHAPTER - 1 


HTTROnJCTI® 


The importance of using porous materials as bearing components 
has been realised because of the fact that the porosity increases the 
stiffness of the bearing It also makes the bearing more stable in 
the dynamic situation The various bearing configurations of porous 
metal bearings have been studied by Cameron, Morgan and Stainsby 
( 1962 ), Joseph and Tao (l966) , Shir and Joseph (1966), H (l972) 
by using both no-slip and slip boundaiy conditions at the porous 
surface The behavior of non-¥ewtonian fluids m a porous matrix 
of a bearing has also been recently studied by Isa (l975) lo. most 
of these studies, it has been pointed out that the load bearing 
capacity decreases with the permeability of the material and with 
the slip coefficients But use of the permeable material increases 
the stiffness 

It has been an usual practice to assume that there is no slip 
at the surface of a permeable material, when there is a flow over it 
But a recent experimental jnvestigation of Beavers and Joseph (l967) 
shows that such an assumption Is at best a crude approximation, to 
the realily This is because of the observed migration of fluid 
tangent to the boundary within the porous material In order to 
account for this Beavers and Joseph proposed a model wherein there 
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is slip at the surface postulate that the slip velocity 

at the pemeable naterface differs from the mean filter velocily 
within the permeable material and that shear effects are trans- 
mitted into the body of the material through a boundary layer 
region hifact the 'adhoc' condition of Beavers and Joseph 
implies that tne slip velocity at the interface is proportional 
to the shear at the wall They argue that the proportionality 
factor depends only on the structure of the material within the 
porous matrix and not on the flow Subsequent eaqjeriments of 
Sparrow, Beavers, Chen and Ilcyd (l972), Beavers, Sparrow and 
Magauson (l970) , Taylor (l97l) for liquids in Poiseuille flows 
and Beavers, Sparrow and Masha (l974) for gas flows substantiate 
the claim of Beavers and Joseph 

The above mentioned results are for developed flows A 
nattaral question would arise as to what should be the condition 
at the surface of a porous ra-erial for a developing flow, as for 
example the entry flow in a porous walled channel One anticipates 
an explicit dependence on a length scale in the direction of the 
flow, the distance from the leading edge playing this role for the 
example mentioned above This is because one expects the migration 
of fluid tangent to the boundary to vary from the leading edge as 
the outer flow will drag more and more fluid inside the matrix of the 
porous blook as the flow advances m the downstream direction The 
mam body of this thesis addresses itself to this question A mathematioal 
nature of the solutions of the problem discussed than would result in 
an explicit formulation of the oonditicxi 



CHAPTER - 2 


A REVIEW OR THE LlTERiTORS 

2 1 Peimeable walls with no slip 

How of an incompressible viscous fluid over a porous material 

has been tne subject matter of several research workers !Ehe 

importance of flows with permeable boundaries seems to have originated, 

specially for the aerodynamical purposes, by ludwig Prandtl (l904) 

An important consideration of a flow over an aerofoil is to determine 

the point of seperation of the flow as the resulting drag thereon 

IS considerable In order to be able to reduce the drag at high 

speeds, it then becomes necessary to deploy some methods which could 

del^ separation Prandtl was able to show that this could be 

achieved to a fairly good extent by using slots in the body and by 

'sucking in' the retarted fluid in the boundary leyer thus ma kin g 

the boundary layer stick close to the body to quite a distance towards 

the trailing edge A natural outcome of such a study was to 

introduce a permeable wall instead of having slots in the body 

Braslow, Burrows, Tetervin and Visconti (1951) have given experimental 

and theoretical results for a flow over an aerofoil with continuous 

7 

suction for a considerably high Reynolds number of the order of 10 
{They use a porous aerofoil to achieve the uniform suction 
Schlichtang(l955) and Rosenhead (l965) give a review of the work done 


in this direction 
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It is to be observed at this stage that laost of vre models 
dealing with tne flow of c. fluid over a porous body xake for the 
boundary ocndition at the wall no-slip and prescribe a suction or 
injection at the v/all But when one coneiders the rarefied flow, 
one will have to give allowance for the slip at the wall 

The effect of wall porosity on the velocity and the pressure 
distribution in a porous walled channel has been studiel by Berman 
( 1953 ) Using a perturbation procedure he is able to deduce that 
the velocity profiles are flatter at the center of the channel and 
steeper at the walls as compared to the Poiseuille flow Yuan (1956) 
gives a further treatment of the above problem lor fairly moderate 
suction Reynolus number Horton and Yuan (l964) have s'^lved the 
entry length problem using Karman-Pohlhausen method Their results 
agree fairly well with those of Schlichting (l955) and Bodoia and 
Osterle (l96l) Recently Raithby and Knuv^sen (l974) have studied 
the i^drodynamic development with suction and blowing The 
governing equation is the vorticity equation and tnqy solve this 
numerically They show that for the case of strong suction, the 
hydrodynamically developed solutions for the velocity profiles 
previously published are not normally attained Their experimental 
results confirm this aspect 

2 2 Permeable Walls with Slip 

Beavers and Joseph ( 1967 ) have givsn some experimental results 
for the flow of an incompressible viscous fluid m a channel one of 
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whose walls is peimeable Because of the- flow of the fluid over 
the permeable material, there is a migration of fluid tangent to 
the boundaiy within the porous body which would not entail one to 
use the familiar no-slip condition at the wall the porous 

wall one means the smooth geometric surface which contains the 
outeimost perimeters of the surface pores of the porous material 
Beavers and Joseph propose a model viiich takes into consideration 
the migration of fluid inside the matrix of the porous material 
by stating that the slip velocity at the boundaiy is proportional 
to the shear at the boundaiy, the proportionality factor depending 
only on the geometric structure and the permeability of the material 
Taylor (l97l) gives a justification for the above model with the 
help of a simple mathematical model and shows that the slip 
coefficient is indeed independent of the geometiy of the measuring 
device Richardson (1971) in a companion paper to the one of 
Tgylor, gives further analysis of the model Saffman (l97l) has 
given a theoretical foundation for the model of Beavers and Joseph 
using a statistical approach He first deduces the Dircy*9 law 
for a non-homo geneous porous medium and using an asymptotic analysis 
he deduces the boundaiy condition of Beavers and Joseph as a 
consequence. 

Sparrow, Beavers, and Hung (l97l) consider channel and Tube 
flows with surface mass transfer using the bomdajy condition of 
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Beavers and Joseph at the wall The analysis takes into conside- 
ration, inj ection/suction at the surface The results show that 
the streamwise pressure gradient may increase or decrease with slip 
as compared to the case when the wall is impermeable Beavers, 
Sparrow, and Magnuson (l970) have given some experimental results 
for parallel flows in a porous walled channel Their results show 
that the mass flow increases where as fricticm factor decreases as 
compared to the solid walled channel flows when slip at the porous 
wall is taken into consideration Sparrow, Beavers, Ohen and Iloyd 
( 1973 ) consider the breakdown of the laminar flow in peimeable walled 
ducts They have conducted experiments with slip velocities at the 
walls inaddition to the case of zero slip The instability Reynolds 
number /for both the theoiy and experiment lie below that corresponding 
to the 3mpermeable walled case 

The effect of velocity slip on the scj,ueeze film between 
rectangular plates has been studied by Wu. (l972) His analysis 
shows that the existence of slip velocity will further reduce the 
load carrying capacity and the response time of the squeeze film 

A recent investigation of Beale and Bader (l974) on coupled 
channel flows with bounding porous walls giv® a rigorous physical 
and mathematical basis for the Beavers and Joseph model They show 
that the slip coefficient a is equivalent to ^'^\i where v is the 
effective viscosity and p the viscosity of the fluid Th^ show 
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"bhai; "the Brmtaiian*s extensioa to Darcy’s law (which takss into 
account a viscous teim in the goveining equations ) accounts for 
the existence of a boundaiy layer region within the po 3 X)US medium 
and that this would be an important factor when one considers 
thin channels Infact, they point out that one need not take a 
recourse to emperical boundaiy conditions when one uses the 
Brinkman’s extension to Darcy's law 

2 5 Plan of the Thesis 

Chapters III, IV and V address themselves to answering aspects 
of flows over a porous flat plate, the entry flow in a porous walled 
channel, and the problem of setting up of convection currents m a 
thin film of fluid over a porous block saturated with that fluid 
respectively For all these problems the condition at the wall 
generally corresponds to the slip condition of Beavers and Joseph 
with a variant The flows considered are all laminar, inccanpressible 
and the fluid is viscous 

It is preferable to use the variant of Beavers and Joseph 
model to the one presented by Neale and Nader (l9T4) for two reasons 
Firstly thou^ the Brinkman equation has received encouraging 
theoretical verifications from several authors, no experimental 
verification seems to have been forthcoming so far Secondly it 
is possibly more appropriate to study the flow in the porous medium 
statistically and match at the boundaiy as definition of point 
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fvinctiona xnsxde the porous medi-um xtself are not rigorous 
Possibly the best course of study would be that of Beavers and 
Joseph Their condition essentially shows the interaction of 
flow m the porous material and the external flow and it could 
perhaps be well-done using a boundary condition rather than 
solving the two flows separately and matching them. 

In chapter III, the flow of an incompressible viscous 
fluid over a porous infinite flat plate saturated with the working 
fluid IS studied The porous plate is assumed to be both homogeneous 
and isotropic As the flow is over a porous substance, as observed 
by Beavers and Joseph (l967), allowance for slip at the surface of 
the material has to be made of This is due to the fact that the 
outer flow drags with it the fluid m the porous matrix near the 
surface As the flow starts from the leading edge, more and more 
of the fluid in the matrix of the porous medium is affecT/ed so liiat 
one expects the slip at the wall to change continuously with the 
distance along the flow direction As in the case of flow over an 
impermeable plate, one looks for similar solutions and as a consequence 
one finds that the slip velocity varies with the square root of the 
distance from the leading edge In addition to the already postulated 
dependence on the shear at the wall The flow equations have been 
solved with this condition as the basis and the results show that 
the full velocity is attained faster as contrasted with the impeimeahle 
plate case 
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In chapter IV, the entiy flow of an incomprebsiDle vxsoous 
fluid in a porous vailed channel is studied Ihe previoua problem 
becomes a motivation for this case in this sense whether the 
groi,/th of the thickness of boundaiy l^er increases or decreases 
as compared with the case of impermeable walled channel One of 
the boundaiy conditions used at the vail is similar to the one 
tnat IS used m the previous problem The method of solution 
consists of using the von-Karman-Pohlhausen 's approach The 
entiy length determined for this problem is shown ro decrease with 
an increasing value of a parameter which is proportional to the 
permeabilily of the porous medium used 

In chapter V, the Benard problem over a porous block is 
studied The model consists of a porous block satur<^ted with the 
woikmg fluid on top of wnich there is a than layer of this fluid 
The bottom of this system is heated so that a dij.ferential 
temperature is maintained across the system Because of the 
interactiaa between the gravity forces and the viscous forces an 
instability sets up m the originally static set up and the aim of 
this investigation is to determine when this happens Tne results 
show that the critical Rayleigh numbers are considerably reduced 
as compared with the case of an impermeable block 



CHAPTER - 3 


LAMMAR PLOW OVER A POROUS PLAT PLATE 


3 1 Introduction 

The problem of the lammar flow of an incompressible viscous 
fluid over an impermeable flat plate was first studied by Blasius 
(l908) Using a non-dimensional stream function he reduced ‘che 
boundary layer equations to an ordinary non-linear differential 
equation and gave its solution using a series method Subsequently, 
various methods have been devised by several workers to get a 
fuller picture of the flow over a flat plate and their results are 
to be found in Schliohting (4§^) and Hosennead (1963) 

One of the subsequent developments in this direction has 
been to study the flow of a fluid over a porous flat plate As 
the plate is porous, fluid can be injected into the external flow 
or sucked from it The chief reason for studying such a problem, 

especially the problan with suction, is that suction causes 
separation to be delayed and this intum will have a stabilizing 
effect on the laminar flow As such it has been customaay to 
associate injection or suction with porous materials prescribing 
these at the wall, the usual no-slip condition at the wall being 


assumed to hold 
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Recently Beavers and Joseph (l967) have proposed an 
emperical relation for the velocity of a flvad at the boundaiy 
for a flow over a permeable body !!liey consider the recti- 
linear flow of a VISCOUS fluid thaxni^ a two duaenaional parallel 
channel formed by an impermeable upper wall and a permeable lower 
wall. The flow is generated over and Inside the porous body by 
the application of a pressure gradient maintained m the longi- 
tudinal darection Because of the flow over the permeable wall, 
there is a migration of fluid tangent to the boundaiy within the 
porous body But it has been oustcmary, prior to the work of 
Beavers and Joseph (1967), to disregard this aspect and use the 
no-slip condition at the surface The no-slip condition m such 
a situation would be the following the tangential velocity is 
zero at the wall whereas the noimial component is prescribed 
Beaveps and Joseph, therefore, postulate that the slip velociiy 
at the interface differs from the mean filter velocity within the 
porous material and that shear effects are transmitted into the 
body of the material through a boundaiy layer region They 
assume that the slip velocity for the free fluid is proportional 
to the shear rate at the boundajy, the proportionality constant 
depending on the geometric structure and the permeability of the 
porous body 331 ^ have tested their model and the experimaatal 


results are satisfactory 
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Taylor (l97l), Hichardson (l97l) and Saffnan (l97l) have 
given theoretical justifications of the Beavers and Joseph nodel 
Taylor gives a simple mathematical model based on the "paint 
brush model" of the above problem He justifies the claim by 
Beavers and Joseph that the slip coefficiait need only depend on 
the structure of the porous material, out not on the features 
of the geometry of the measuring device Richardson in a 
companion paper gives further analysis of the Taylor model 
Saffman gives a justification of the proposed condition employing 
a statistical approach He first deduces the Darcy's law for a 
non-homogeneous porous medium and using an asymptotic analysis 
he deduces the boundary condition as a consequence 

In the present investigation, the flow of an incompressible 
VISCOUS fluid over an infinite porous flat plate is restudied 
As the plate is porous, the no-slip condition is discarded 
Instead, a boundary condition similar to that of Beavers and 
Joseph IS assumed The method of solution is the numerical 

solution 

3 2 The Model 

The physical problem is the lammar flow of an Incompressible 
viscous fluid over a flat porous plate The plate is assumed to 
be both homogeneous and isotropic The plate is saturated wxth 
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the fluid before the flow starts Ihe x— coordinate is tahen 
along the plate and the y— coordinate normal to it with the 
leading edge of the plate being taken as the origin A sketch 
of the model is as shown m 5h.g 3 1 There is a uniform 
flow of the fluid with a velocily upstream of the plate 
The problem is to know how the velocity developes on the plate 

The governing differential equations are 


u 


8x 


+ 




3y 


(3 l) 


8u 

8x 





(5 2) 


where v is the kinematic viscosity of the fluid The 
boundaiy conditions at the surface are 

(|)( — — ) ^ = — u,v = 0 at y = 0 (3 3) 

u as y « ( 34 ) 


Here a is the slip coefficient which depends on the geometry 
of the material and k is its permeability k is a measure of 
the fluid conductivity of the material 

Some comments about the boundaiy conditions as to their 
forms are relevant here The first condition corresponds to the 
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slip condition of Beavers and Joseph, £ind as li their paper 
IB taken to be an emperical relation ?/hen the full velocity 
profile IS attained, which in practice is attained for a finite 
distance x, the resulting condition at this station and 
onwards reduces to the condition of Beavers and Joseph Ihe 
adhoc condition assumes that the slip is a function of x, 
the distanf’u from the leading edge, till the full profile is 
attained As there is no preferred length for the model, 

' f 

— / has been taken as a length scale 

V/ 


As in the ease of BlasiQ.'S' problem, when one searches for 
similar profiles for the flow, the function ij) becomes 

/ — — Ihe second condition m Equation (3 3) shows that there 

IS no suction or inj ection fhis condition differs from the 
usual analyses of flows over a porous uoJy wherein v is 
prescribed 


Introducing a non-dimensional stream-function f(h) and 
the transformations 


=5 yy — 


yv — , ipr= f(n) 


Equations (3 l) aiid (3 2) reduce to 

orilA. „lfS ^ it 
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2 f" + ff" =0 (3 5) 

where primes denote differentiation with respect to p 

As similar solutions are sought for, the boundaiy 
conditions become 

f(o) = 0, f"(o) = ef'(o) , f'^ 1 as n (3 6) 

where use has been made that the function <J> can be only 

/ xUoo yoj 

/ __ 0 _ 

3 3 Solution 

The method of solution is the numerical integration IThe 
differential equation (35) has been integrated using the Eunge- 
Kutta method (Oamahan, Luther and WiUces (l969)) with the 
relevant boundaiy conditions An alternate method for solving 
the Equation (3 5) is also given 

The differential Equation (35) can be written as a set 
of three ordinaiy differential equations 


dn “ % 

_ 1 
db ~ " 2 


(3 7) 
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The boundd,3y conditions become 

g^(0) = 0, g^(o) = gg^Co), gg -*• 1 as n -o (3 8) 

The above system can be treated as 'mitial' value problans 
for a given initial condition g2 (o)j as tv/o of the initial 
conditions could then be had Then a search for a value of 
g^ at b = 0 IS made that will generate a solution for which 
g^ 1 as n “ This can be accomplished by having an 

iterative procedure for g2(8) Then the whole system of 
equations can be solved using the Eunge-Kutta method 

The range of values of 3 is from 0 to <» , the former 
value corresponding to the infinite permeability case-hence an 
infinite slip- which is an 'ideal' situation and "che latter 
case corresponds the Blasiue problem Thus the velocity curves 
are to lie between the Blasius profile and the 'ideal' profile 
f = 1 

The velocity profiles for some representative values of 
e are drawn and are shown in Pig 3 2 

An alternative approach for solving the Equations (35) 
and (3 6) can be given The advantage of this method is that 
the given boundary value problem is reduced to another boundary 
value problem whereon a search for the initial condition need 


not be made 
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Making the transformations 

r)* = and f = 2 Bg (3 9) 

the given boundaiy value problem reduces to 

g"’ + gg" ==0 (3 10) 

With the boundaiy conditions 

g(o) = 0, g"(o) = g'(o), g' n* “ (311) 

where primes now denote differentiation with respect to p* 

The above problem could be solved as in the previous 
problem with this difference The integration can be done 
starting with some initial condition, the integration being done 
till g' attains a constant value which will give the value 
of B A sketch of the profiles g' against ri* is given 

in ihg 5 3 

3 4 Conclusions 

The velocily p 3 ?of lies for different values of the 
parameter B are shown in Pig 3 2 The Blasius profile 
corresponds to the value B 'equal' "to infinity It is to be 
observed that the profiles for increasing values of B are 
one below the other Hie curves for smaller values of B have 
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sharper curvatures and on already observed, the ideal situation 
of infinite permeability corresponds to the curve f = 1 
Ihis suggests that the full velocity is attained faster as 
compared to the Blasius profile 

It IS to be noted that the usual compatibili-ty conditions 
are not available for the above model because of the allowance 
for the slip at the wall and the absence of any suction or 
inj ection at the wall As such no quantitative comparisons 
can be made as the problems are physically different But 
some remaiks could be made about the similarities between the 
model considered and the model describing the flow over a 
porous flat plate with suction vathout slip The asymptotic 
suction profile is sharper as contrasted with the Blasius 
profile ( Schlichting, pp 270 ) ihr the problem under 
consideration a similar result is true The measuremaits of 
Head (l95l) seems to suggest similar results as he has made 
measurements for relatively small suction and his results being 
true from a short distances away from the leading edge 
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CHAPTER - 4 


ETTRI PEOW IH A POROUS WAILED OxiADE'JEL 


4 1 Introduction 

The flow of an incompressible viscous fluid in a channel 
con generally be divided into two regions an entrance region 
which is characterized by the changing cf velocity profiles 
continuously in the direction of the flow and a fully aeveloped 
region wherein velocity profiles exist independent of the flow 
direction 

The steady laminar flow of an incompressible Hewtonian fluid 
in the entrmce region or the inlet section of a channel nas 
been investigated by Schlichting (l955), Bodoia and Osterle (l96l), 
Collins, Morton and Schowdlter (l962) ana others These authors 
use the usual boundary layer theory In all these cases the 
velocity distribution has been assumed to be flat at the entrance 
of the channel and the boundary layer approximation of neglecting 

the viscous term y •~:r and the pressure gradient nomal to the 

9x^ 


flow direction has been made Scnlichting (l955) using the original 
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Blasius (19O8) technique for a flat plate obtained the velocity 
profile for the entrance region by means of a roatching procedure 
between an upstream and a downstream solution to the differential 
equation governing the system Bodoia and Osterle (l96l) solved 
the entiy problem by means of a finite difference procedure and 
obtained results which differed firom Schlichting's solution 
They explained the discrepancy as follows firstly there is 
always a discontinuity introduced at the point where the upstream 
and downstream solutions were matched and secondly Schlichting’s 
assumption of the second derivative of the foiward velocity with 
respect to the cross channel position variable being zuro is incorrect 
Oollins and Schowalter (1962) used Schlichting’s method with 
refinements Eecently van I^ke (l970) and S D R Wilson (l97l) 

have worked on this problem of development of c. channel flow using 
essentially the same techniques 

The effect of wall porosity on the velocity and the pressure 
distribution in a porous walled channel has been investigated possibly 
for the first time by Berman (l955) As ne studies the fully 
developed two-dimensional flow, he uses a stream function to reduce 
the full Favier-Stokes equations to a third order non-linear 
differential equation with the usual boundary conditions with a 
prescribed constant suction at the wall Using a perturbation 
procedure he is able to deduce that the velocity profiles are 
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flatljer at the centre of the channel and steeper at the walls as 
ccanpared to the Poiseuille flow and that the pressure drop in 
the direction of the flow is found to be appreciably less in the 
porous walled channel corresponding to an impermeable walled 
channel of the same dimensions and the same entry Ra 5 molds number 
Yuan ( 1956 ) gives a further treatment of the above problem when 
the suction Reynolds number is fairly moderrate Horton and 
Yuan (1964) have solved the entry flow problem using the Karrn^- 
Pohlhansen method Ihey give solutions for both similar as well 
as non-similar profiles in the entry region The parameter used 
m the calculation is a Reynolds number depending on the suction 
velocity at the wall Their results agree fairly well with those 
of Schlxchting (l955) and Bodoia and Osterle (l96l) 

Sparrow, Beavers and Hong (l97l) have studied the problem 
of channel and tube flows with surface mass transfer and velocity 
slip at wall They use the Beavers and Joseph model for thu 
boundary condition at the wall Their analysis shows some 
interesting results because of fhe presence of slip at the wall 
They show that relative to tne situation viiere there is no slip, 
the streamwise pressure gradient may either increase or decrease 
in the presence of slip and m^ even undergo a change of sigi 

However, the entry flow in channels with permeable walls 
allowing slip does not eean to have been studied 30 far 
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In the present Investigation, the fntiy fluw of an 
mcjmpressible visgcus fluid in a pemeable walled choinsl is 
stuaied The peimeable walls are assnaea to be hon^geneous 
and isotroi lo and they are saturated with tne workin g fluid to 
start with Generally, two types cf flovi at the entiy re 
assuned (i) the flow is unifoim at the entra.ice of the channel 
and (ii) the flov/ upstream of the leading edge of one of tne walls 
is assumed to be uniform Assumption (i) is used more frequently 
though the second assumption is more realistic piiysically Wang 
and longwell ( 1964 ) have shown that very definite effects are 
transmitted upstream from the leading edge of the plates at a 
Reynolds number of 300 for the first tyj e of initial condition 
They have also found that the overall pressure drop and the inlet 
length predicted by boundary l^er theory methods are in reasonable 
agreement with those of the boundary conditions corresponding 
to the case (ii) They feel that this night be due to compensating 
errors caused by the several assumptions made m the boundary 
layer theory von lyke, M (l970) h'^s observed that use of condition 
(i) introduces a weak virticily at the inlet and that this leads 
to certain complications 

In the present analysis the flow is assumed to be uniform 
at the entry The method of sulution adopted is based on the 
van Kam^-Pohlhausen method Briefly the method is as follows 
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The boundaiy layer equations are integrated once to obtain the 
Karm^ momentum integral equation A form of velocity profile 
u(x,y) IS then sought which satisfies the momentum equation and 
the relevant boundazy conditions «/hen tie assumed velocity 
profile IB substituted on the momentum integral equation, these 
results a differential equation for the shape factor v/nicn is a 
function of the external pressure gradient and the boundaiy layer 
thickness solving for the boundaiy layer tnickness, the 

other variable like the displacement thickness, the momentum 
thickness are determined Pohlhausen used a quartic expansion 
for the velocity profile 

Several improvements of the Karman— Pohlhausen method are 
available Mangier (1944) suggested the use of a polynoiaial of 
the n—th degree for the velocity profile which automatically 
satisfies some of the boundaiy conditions Tinman (l949) has 
suggested the use of certain indefinite integrals for the velocity 
profile Several other 'approximate’ methods are available which 
use m addition to the momentum integral equation, the energy 
Integral equation also All these methods work fairly well when 
the flows are accelerating - the flow developing in a channel is 
an example of such a flow A fairly exhaustive treatment of 

these approximate methods is available in Rosenhead (1963) and 
Curie (1962) 
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4 2 Model and Mathematical Ihrmulation 


The model consists of a two dimeasioncil steady ilow of an 
incompressible viscous fluid in a peime-ble wallea charnel The 
Walls are 2h units apart cJid ^re assiuaed to oe botn homogeneous 
and isotropic The walls of the cnannel are asburntd ro be 
satur<ated with the working fluid to begin \/ith Tne permeaoilitj, 
of the channel walls is k Because of the symmetry in the flo ; 
field, only the flow in one half of the channel is studied 
A coordinate frame is fixed at the leading edge, it being taken as 
the origin The x— axis is along the flow direction along the 
cnannel and the y-axis is m a transverse direction(see Big 5 1) 


The governing equations of motion are 


u 


3u 

9x 


+ V 


3u 

oy 


ax 


(p/p) + '^ ( + 

3X 



(4 l) 
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) 
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Jli + II « 0 (^5) 

ax ay 


The boundary conditions are 




f(x) 


3y 


V = 0 


(4 4) 


at y = 0 


u 


a 
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y = 5 


lii 

3y 


2 

a u 


= 0 , u = f 


(4 5) 


Hore U xs the central line velocity, k. is the permeability 
of the porous walls, a is the slip coefficient and 5= 6(x) is 
the local boundary layer thicknoss 


Let 


u = 


1^ 

U 


vh 
U (5 


2 . 

■«! 


-,5 = ^ 
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Pe =5 — ;j- , where is the constant external velocity, 2h is the 

width of the channel 


Then maring the Prandtl approximations, these equations 
(4 i), (4 2 ) and (4 3 ) become (after dropping bars) 


^ SJi y Bii _ 

ox On 



-L 

Re 




(4 6) 


Ou 
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= 0 


The non-dimensional boundary conditions are 

f(x) 2_u 


at ri = 0 


u 


On 


, 'U'= 0 


n = 1 


Tu, 

9n 


o^u 

9n^ 


0 , u = U 


(4 7) 


(4 8) 

(4 9 ) 


Integration of the equation (56) from n = 0 to n = 1 ond use 
of continuity equation gives 
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1 5uj 

Re ^2 t)Ti [n =0 



u(u-u)dn 



1 

I 


o 




(4 10) 


11; IS 1;he well Itiiowel Kanuan^s Integral equation 
4 5 Solution 

N 

The method of solution adopted is similar to the von-KainiM— 
Pohlhausen method The velocity in the boundaiy layer is expcuaded 
as a quartic in the non— dimensional transverse variable n The 
coefficients are determined by using sufficient number of boundary 
conaitions from equations (49) let then the velocily profiles 
be giren by 

■^ = I t <f>(x) I , (4 11) 

0 0 

where <l)(x) is s function of x alone and a. and B are constants 

1 1 

The constants in the above expansion depend upon the degree of the 
polynomial chosen The form chosen as in Equatiob (4 11 ) will be 
useful later when discussion for the similar andnon-similar profiles 
IS done As for example, in the case of similar profiles has 

to be a function of n alone so that the second term in the right 
hand side of Equation (4 1l) is dropped for that discussion 
Discussion for the quartic polynomial alone is done in the subsequent 


sections 
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Solution for similar profiles 

If the assumption is made tnat the velocity profiles in 
the boundary laj er are similar at all channel cross-sections, 

u 4 ^ 

then Equation (4 1 1 ) will reduce to ^ ^ An This ctlso 

u 0 ^ 

f (x ) / \ 

implies in turn that the coefficient “ ^ Equation (4 8; 

be a constant say, d The coefficients A^’s have to satisfy the 
Equations (4 8) and (4 9) when 

A^ + Ai + A^ + A3 + A^ = 1 

A^ + 2A^+ 3A3+ 4A^= 0 
2 A 2 + 6A3+12A^= 0 
6A3+12A^= 0 
"o = 

Thus 

. 4d , _ 4._ ^ §— 

\ " 1 + 4d ' ^1 " 1 + 4d » 2 1 + 4d 

h ‘ “ -TTw 


Hence 


u 

TJ 


1 

(1 + 4d) 


[(1 + 4d) - (1 



] 


(4 i 2 ) 
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ilso 

7 


1 

/ u dn = 
o 


20d + 4 
5(l + 4d) 


1 

/ 


0 


2 

u dn 


(720d^ + 288d 4-32) ? 

— i 

45 (l + 4d)^ 


(a- 13) 


(4 14) 


The central line velocity U which is a function of x 
alone needs to be deteiminai Tnis Ccin be done since n definite 
form of expression for u has already been assumed satisfying 
the boundaiy conditions Also at two different sections of the 
channel^ the inflov/ across one section must be equal to the 
outflow across the other section of the channel So much so, the 
inflow across the section at the entrance must be equal to the 
outflow across any section of the channel (The flow inside the 
channel itself could be divided into two halves-flow m the boimdary 
layer where thi^ velocity component is u and tne flow outside the 
boundary layer, the accelerating core, where U is the velocity 

Thus writing the mass conservation in a dimensional form 
one gets 

2h 2h 

/ u dy = / Ucj dy 
0 0 


(4 15) 
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or non- 


dimensionally C -btst-v'-i ^ 

1 h/& 

j udh+s/ Udn=h 
0 1 


(4 1 6) 


3tituting ifi Equation (4 16) for j u dn from Equation 


(4 13 ) and simplifying one obtains 


1 


where 


(4 17) 


B = - 


5(1 + 4d) 


Ihus the momentum integral equation becomes 


^ ^2 ^ 2(20d+ 4) (720d + 288d + 32 ) 

<52 (i+4d) 5(1 + 4d) 45(1 +4d)^ 


( 20 d + 4 ) 

+ 1 ] 

5(1 + 4d) 


(I08d + 17) 

45(1 + 4d)^ 


i- A ^ 

He g2 dx ’ 


(4 18) 


where A = 


I08d + 17 
45(1 + 4d) 
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Substitutxag for — from Equation (4 17) and integrating 
from 0 to jc, 

2 

x=-_ [__^_2Jlnul (4 19) 

4B ^ 

Working on similar lines for the case of a rigid boundary j 
that is, #ien u = 0 at n = 0, one gets 

2 

X = ^ Ee - - 2 «.nu ] (4 20) 

This equation turns out to be a particular case of Equation (4 19) 
for d = 0 

4 4 Eon-Similar Profiles 

For the velocity profiles in the boundary layer to be 
non-similar, the velocily profile is as given an Equation 
(4 11) As only quartic polynomials are being considered 

I + 4 >(x) I (423) 

0 0 

Without loss of generality, can be taken to be 1 Some 

of the constants A^, B^ can be found as in the previous case 
Thus using the Equations (4 8) and (4 9) one has 

■^0 ^ \ ^ ^ ^“^4 

®o “ ■“ '^®4 " ^ i \ ~ 684. 
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Equation (48) yields 


or 


Aq + = g(A^ + <|)Bp 

, _ 1 

1+ 4g ’ 


(4 2^) 


where g 


^ t(^) 

a 6 


» A^ + 


(4 25) 


!rhe function g(x) treated as a function of the parameter < 
will have tne range (O, ) , the end points corresponding to the 

rigid and the 'completely' permeable case (Ehe latter case 
IB ptysicalJy non -plausible Thus the range of the function ilJ 
IS [ -1, O) 

Equation (4 25 ) then reduces to 


u 

U 


(l +4^) - 4 4^n+ 6x[it)^ ~ 44'n^ + 


= 1 + 4' ( 1 - n) 


4 


(4 26) 


Jilso 


1 

/ u dn = II (l + ^ ), 


1 

/ 


o 


U^dT) 





(4 27) 


(4 28) 


Erom the boundary condition given by Equation (4 8) one has 
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4 4 at n = o 

3x dx Re ^ 2 

c 6r) 


(4 29) 


As in the previous case U is determned by -che use of the 
conservation law 


1 ,Vs 

5 / udh+d/ Udn=sh 

o 1 


(4 30) 


Subs-cituting lor the former integral 


U - 


(S 

^"5 E 


(4 31 ) 


Equation (4 29) then gives 


J2 h = 1 + 

Re 52 dx ^ ^ dx dx ^ 


(4 32 ) 


Substituting for the integrals m the moma turn Equation (4 IO), 
one further obtains 




Re 52 


5 9 


541 dx 9 dx 


(4 33 ) 


ELiminating —r between the above two equations 
Re §2 


dU _ dx 
d^; ^ 

dx 


U / 6 -I- 5 t|> \ 

"4; ^ 3 + 54' ' 


(4 34 ) 


On integration, 
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(5^ + 3) 


U = c 


(4 35) 


r 


where c is a constant of integration 


Equatioi (4 32) with the help of Equation (4 3l) and (4 35) 
would then become 


dx _ 23 Re 

dll' “ 12 


(9 + 8’!') c (5’f' + 3)]^ 
(3 + 5 ’f’ )^ 


(4 36^ 


The above equation can be integrated m a straight forward 
fashion Before doing so, certain commenteare due 

The flow in the channel is an accelerating flow so that 


12 

dx 


> 0 But 


22 , - ^ 12 

dx “ di|; dx 


(4 37) 


T TVni fi — and need to have the same signs But from 
d4» dx 

dx 

Equation (4 36), for the range of f which is [-1>0), — > 0 So 


much so, > 0 But 


= ^ ( 5 ^+ 6 ) 


(4 38) 


Hence o is positive 
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Also since U IS positive j and c is dso positive j 
Equc*tion (^t- 35) shows that 


(5 4'+ 3) > 0 
And thus 

- 6 < < 0 (4 39 '' 


Integrating Equation (4^36) one gets 


= _ ^ RefirZOc on ( 54>+ 3 n _ 9 cf 8 

12 9 ^ ^ 4 V 3 ,3 

4 4 * 


4 2 £ _ ^ 


4- 


5(3+54') 


] + A., 


(^ 40) 


where A is a constant of integration 


Since L = 1 at x= 0, from Equation (4 35) one has 


4'q = c (3 + 54'q) j 


(4 41) 


where 4' is the value of 4- at x = 0 
0 

4 5 Discussion 

The entry length for the channel could be defined in two 
ways (i) through Itie boundary layer developing in the channel 
and ( 11 ) throu^ the fully developed velocity profile Considering 


the former case 
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Deiarition I 

(Jhe entiy length for a porous \/all2d cdcr' lel is th-t 
length along tne channel when the bouadaiy laver grov/s asympto- 
tically to half the channel width 

5 y sj'tmaetxy of the flow in the cnannel vath respect to 
the center of the chcinnel, the boundaiy layers developing on the 
two walls of the channel meet and the distance of this point 
from the entrance is the required entry length This sort of 
definition would give a lower boxind for the entry length, (see 
I) vies (1972)) 

But a more appropriate definition would he 
Definition II 

The entry length for a porous walled channel is that lengtn 
along the channel attained when the veiocily profile of the 
developing flow m the channel attains the value of the velocidy 
of the fully developed flow t the center of the chainel 

The procedure for the determination of the entry length 
for the similar profiles in the boundary l^er is c*s follows 
when the profiles are similarj the boundary condition given by 
Equation (4 8 ) is to reduce to 
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u 



n = 0 


(4 42) 


But for the fully developed c 3 St, 1 e , wacn the dineiisional form 
of Beavers and Joseph boundary condition is used 

(4 43 

or in the non-diraensional form 
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3u 
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at y = 0 
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E ii ^ 

ah 6 9p 


at n = 0 


(4 4-4 ) 


If Xx is the entry leng-th, then 
u 


a = £ „ K ' K= 

ah 








(4 45) 


But then from Equation (4 17 )s 
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5(l+4d) 



5(l+4d) 

5(l+4d) 


(4 46) 


Using the boundary condition given by Equation (4 43 )> 

show that u , the fully developed value of tne c^-ntre line 
max ' 

velocity is 


max 2 M + 3K ^ 


(4 47) 


Thus using the Definition II of the entry length, one has 
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d = 


- 5 + *^25 + 240K (1 + ~2K) 
40 


(4 ^8) 


The entry length for the Definition I is calculated hy 

taking 5/h = i m Equation (4 46) \/hich determines D Tnen 

X could be determined from Equation (4 19) for different 
L 

values of d 

The entry lengths corresponding to different values of K 
tevebeen calculated for the above two cases and they are tabulated 
in Table (4 1 ) 


Eor the non-similar profiles, a similar approach would give one 



-5 + *^25 + 240K (1 + 2K) 
40 


(4 49 ) 


Butt 



1 

1 + 4g(x^) 


(4 50 ) 


using the Definition II of the entry length, 


3+5 


/ 1 + 2 ^^ \ 
^ 1 + 


(4 51 ) 


Given a K, 'I' (x^^) coxad be determined from Equation (4 50 ) and 
hence c from Equation (4 5l) "^alue of 4' (x) at 

X = 0 could then be determined from Equation (4 41 ) 
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The entiy length corresponding to the given K then 

ho 

could be determined from Equation (4 

The entry length corresponding to Definition I, i e , in 
terms of the boundaiy layer growing to half the channel width is 
calculated as in the previous case of similar profiles j since one 
has a relation between U and ^ /h. from Equation (4 31 ), a 
relation between U aid f from Equation (4 35) and a relation 
between, c and in Equation ! (4 41 ) 

A tabulation of entiy lengths x_ for different values of 

Xi 

K are given in Table ( 42 ) for the case of both the definition 

For both the similar and non -similar velocity profiles in 
the boundaiy layer of the developing flow, one notices from the 
Tables (4 I) and (4 2), that the entiy length reduces wiHa the 
increasing values of the parameter K, 1 e , with the permeability 
of the walls for a given value of a , the slip coefficioat The 
chosen values of K are relatively small compared to the range of 
Z vAxLoh can be ( 0, “ ) Also one notices that the value of the 
entiy lengths for the definition given inteims of the boundaiy 
layer attaining half the channel width is less than that of the cri- 
terion In terms of the value of the fully developed profile which 
conforms with the statement of Davies ( 1972) 
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Table 4 1 

Entiy length Similar Profile.^ 


K Eatiy Length 



Defn 

I 


Defn. 

II 




-3 


-2 

0 

8 972 

X 

10 ^ 

5 513 

X 10 




-3 


-2 

0 001 

8 750 

X 

10 ^ 

5 283 

X 10 




-3 

5 277 

^-2 

0 002 

8 734 

X 

10 ^ 

X 10 




-3 



0 003 

•8 718 

X 

10 

5 270 

X 10 




-3 

5 263 

“2 

0 004 

8 712 

X 

10 ^ 

X 10 




-3 

5 257 

“2 

0 005 

8 685 

X 

10 ^ 

X 10 

0 006 

8 669 

X 

10-5 

5 250 

1 

o 

X 


0 007 
0 008 
0 009 
0 010 
0 015 
0 017 
0 018 


8 653 10“^ 

8 637 10"^ 
8 622 X 10”^ 

8 606 X ^o~^ 

8 527 X 10”^ 

8 496 X io'^ 
-3 

8 481 X 10 


5 244 X 10 
5 237 X 10' 
5 230 X io' 
5 224 X io' 
5 190 X io‘ 
5 176 X 10' 
5 170 X io’ 



0 019 


5 163 X 10 
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Table 4 2 

Hatry Leagth Non-Soiailar Profiles 


K 




Datiy length 





Defn 

I 


Defn 

II 

0 35 

8 

435 

X 

io“^ 

5 

872 

-0 

X 10 “ 

0 40 

3 

956 

X 

10-3 

2 

734 

X 10"^ 

0 45 

2 

387 

X 

,0-3 

1 

637 

X 10'^ 

0 50 

1 

615 

X 

10-3 

1 

099 

^“2 

X 10 

0 55 

1 

164 

X 

10-3 

7 

885 

X 10”^ 

0 60 

8 

792 

X 

10-^ 

5 

917 

X 10'^ 

0 65 

6 

850 

X 

10-^ 

4 

588 

X 10’^ 

0 70 

5 

463 

X 

lo"'^ 

3 

647 

“3 

X 10 

0 75 

4 

448 

X 

10~^ 

2 

951 

-3 

X 10 

0 80 

3 

681 

X 

-4 

10^ 

2 

436 

-3 

X 10 

1 00 

1 

931 

X 

-4 

10 

1 

265 

-3 

X 10 

1 10 

1 

473 

X 

10-* 

9 

595 

“4 

X 10 

1 20 

1 

150 

X 

lo""^ 

7 

463 

X , 0 “* 

1 50 

6 

118 

X 

,0-3 

3 

930 

X 10"^ 

2 0 

2 

712 

X 

,0-3 

1 

718 

X lo"-* 
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Entry flow w h rorous walled channel 






chaptur - 5 


BMTiffiD PROBLiM \/IIH A POECuS BODBDiiHr 


5 1 Introduction 

The classicdl Benard problem — the setting up of convection 
currents m a thin film of an inccmpressible viscous fluid l^er 
originally at rest due to an adverse temperature gradient across 
it - has been subject to a thorough study by a number of autnors 
(Chandrasekhar (l96l)) The basic mechanism for the instability 
of the system has been clearly explained by Pearson J R A (1958) 
in terms of surface tension forces while a similar paroblem of the 
breakdown of a fluid subject to a vertical tempercture gradient 
in a homogeneous porous medium and the ensuirg convective flow has 
been studied by lapwood (l948) Lapwood, i sing the linear stability 
theory and the Boussmesq approximation arrives t a fourth order 
differential equation for the temperature perturoation e As 
regards to the boundary conditions at the upper and lower surfaces of 
the porous block, the temperature has to satisfy one condition at 
each surface Bit since the solution of the fourth order differential 
equation admits four degrees of freedom, the velocity can be made to 
obey one condition only at each boundary vdiich implies that allowance 
for a state of no-slip cannot be made as this would require two 
conditions lapwood suggests then that the solution of the differential 
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equation is Yu.lid at best upto the nisighbourhood of the botmdaiy 
and oeases to appfy in a boundaiy layer 

It has been customaiy to assume that the tange^itial relocity 
IS zero at tne surface for a flo# over a porous suiface Beavers 
and Joseph (196?) have recently proposed an adhoo ’slip-conaition’ 
for a flow of the type as mentioned above 

The validity of the model proposed by Beavers and Joseph 
has been tested experimentally and theoretically for certain flow 
configurations 

The linear stability of a laminar flow m a parallel plate 
ohannel, one of whose walls is porous, has been investigated by 
Sparrow, Beavers, Chen and Lloyd (l973) fh^ use the 'adhoc' 
boundary condition of Beavers and Joseph The validily of this 
condition has been established by faylox {l97l)> Richardson (l97l) 
in companion papers and by Safftaan (l97l) Sparrow et al observe 
that the behavior of the critical Reynolds number depends m a 
complex way on the non— dimensional permeability A of the porous 
medium which oan be regarded as a measure of the extent of the 
interaction between the flows in the channel and the porous material 
In certain ranges of A , it is observed that the critical Reynolds 
ntamber is substantially reduced while in otners it increases An 
explanation lo advanced m terms of several mechanisms by which the 
permeable material affects the laminar stability lamt of the channel flow 
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At places where the curve shows a decreasing tendenc/, it is stated 
that the non-vanishing disturbance velocity is believed to dominate 
while m the range in \/hich the curve shows an increasing tendency; 
the effect of thi- sicewness of the main velocity profile lonundtes 

The m.^in purpose of the present investigation is to understand 
this explanation more clearly If the main velocity distribution is 
totally absent, then how the non-vanishing of the disturbance velocities 
is going to affect the stability criterion ’ Wilx it still be 
destabilizing for all the values of the parameter A ’ An answer is 
sought for these questions 

5 2 The Physical Problem and the Basic Equations 

The model consists of a rectangular porous block of thickness 
h whose underside is a rigid plate The porous medium is assumed 
to be both homogeneous and isotropic There Is a thin layer of an 
incompressible viscous fluid of thickness d on top of the po3rous 
material which is saturated with the same fluid %• heating the 
underside of the porous material, an adverse temperature gradient is 
maintained across the system Because of theimal conduction the 
density of the fluid in the porous medium as well as that of the 
fluid layer changes This will lead to an instability whenever the 
adverse temperature gradient exceeds a certain value A sketch of 
the model is an shown m Pig 5 I 
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The govexnang equations oi motion for tne fluid are the 
usual Navier—Stokes equations and the heat conduction eouation 
For the fluid m the porous medium, they are the unsteady form of 
the Darcy’s law and the modified heat conduction equation (Katto 
and Masuoka, 1967 ) In cartesian tensor notation tney arc 

For fluid 
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For fluid in permeable medium 


(5 1) 
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0?he aiaertial terms in the Darcy's law have heen neglected aS they 

are very small compared to the other terms Here £ is the 

porosity ol the material, A* the permeability, k the coefficient 

of thermal —11 c the specific hi^ar at cor start 

P 

pressure The subscripts f and m reier to the fluid medium 
and the porous material-AA-hilvijAAljl 

The equations of state have to be supplemented to this set 
of equations and they are 

p = P [ 1 ~ a(T - T )] 

^ ^ ( 52 ) 

p* = Po f 1 - - Tj] 

where a , are the coefficients of volume expansion of thw 

fluid cjnd of the fluid in the porous material, is the 

temperature at the interface 

5 3 The Perturbation Equations 

The initial state when there are no motions is governed by 




where xs the temperature at z = 0 Oie notes that the 

temperature is continuous at the interface z = h She values of 
0^ and. 6^ are right now unknown But if the temperature T 2 
at z = h + d IS prescribed and the temperature ilux be continuous 
at the interface z = h then one ootains that 


r~ ''T - 

[ K ~ = K - — at z = n 3 

f 9 z m 9z 


T = T - B„h - 6„d 
2 o m f 


0 _ Tc 0 Xfn 

f f mm-' 


(5 't) 

fTier-meJ- CoTiAuc-t 


These equations deteimine 3^ and 3^ as 




_ "a’ 

^ (hKj + d kJ 


(5 5 ) 


1? (t - T ) 

— HI o ^ 


Also the density at the interface is to be continuous so that 


1 e 


D = P [ 1 - a (I,- I) 3 
10 110 


p =sp (i+a, Sh) 


(5 6) 


Under the Boussinesq approximation P., -Pq 
is a fluid of density saturating the peimeable medium and lying 
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over making a layer of thickness d The system now is heated 
from below maantaiumg the top and the bottom w 11 s at temperatures 
Olid (> 12 ) Through conduction the densily gets modified 
The purpose now is to investigate the onsex of convection in the 
fluid layer 

The flow variables are now perturbed in the usual fashion 
The linearized perturbation equations under xhe Boussinesq approxi- 
mation are 


For fluid 
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For peimeable Medium 
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T 116 Isst "tsiiu "fclio first s^uation of tbs Sst (5 7 ) can bs diroppod 
as there occurs the product of a and which are both -■=»»"« i1 

lOie physical condition at the upper ngid plate adjacent to 
the fluid, requires that the velocily be zero However at the lower 
rigid plate adjacent to the porous material (1 e at z = O) only 
normal component of velocity can vanish At the interface z = h, 
transverse velocity, temperature, normal stresses and temperature 
flux are to be csontinuous JUrther to include the interaction 
between the fluid and the porous medium, the conditions proposed by 
Beavers and Joseph are made use of !rhen the boundary conditions 
can be written as 


At 

At 

At 


z = 0 
z = d+h 

z = h 


0 ' 

0 


=s 0 = w' 


= 0 = w = 


3w 

9z 


, - 30' 

W = w',0 = 0 ', 


J 51 

f 3z * 


(5 8 ) 


- p’ = + 

I- (u,v) [(u,v) - (u', v‘)3 

3Z /J* 

Here u,v and u', v' are the velocities of the fluid and the 
velocities of the fluid in the porous medium, a ^ is the slip 
coefficient u* j v* ^ determined from the 3)ar<y*s law 


THow using the normal mode analysis the perturbed flow 


variables ar*e expanded as 
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f(x, 3 f,Zjt) - f(z) exp [ik X + ik y + pf ] (5 9) 

y 

where p is complex in general Making the transformations 
z = T » a =® M, 0 = - , p = , p = 

d * ’ V ’ m K ’ K„ » 


A = 


.,hs|-,k"=k".k2,v=Ji. 

,2 dz ’ X y » p 


(5 10) 


“ '•Q 

the perturbation equations after simplifications and eliminations 
reduce to 


(d^- a^) (d^- a^-o) w = a^ e 

3/ 


(D^- a^- op^)e = - ~ d‘w 

tCj 

n « R K Aa^et 

((7 At 1) (I|2- a®)»' = - -g ' S - Jg 


6 


|< “11 


(5 11) 


f ,2 


(D^- a^- oP^) e' = - 




m 


0< z < I 


w' 


(5 12) 


m 


gaa 




1'«“l “in ^4 
a 


where R_ = — — d , R = — 

f <fV ’ m K^v 


lUrther^ perturbations inpressure within fluid and porous medium 
will be given by the following 


P = (d^- a^*- a ) m 

a'^d 


>»_ 


eAda^ 


a»% 

AoCr No. ^ «!»w 


suss 


(l +0A) DW' 
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Making use of tiie above expressions for p and P', the boundaiy 
conditions (5 8) came out to be 


At 

z 

= 0 

e‘ 

0 

II 

II 


At 

z 

== ■> +1 

e 

H 

II 

0 

11 

(5 13) 

At 


h 

0 



z 

^ d 

= 6', w = w', K D e* = K„ D6 
' m f 



(d^- 2a^-a )Dw + (l + CfA) Dw' = 0 
eA 

2 “o 

D w = (Dw - Dw' ) 

ViTA 

5 4 Evaluation of the Critical Eayleigh numbers 

The system of equations (5 14), (5 15) along with the set 

(5*13) form eigen value problans for the Rayleigh numbers R^ and 

R Hie critical value of these numbers will be determined by the 
m 

use of Galerkin’s method 


In tile following discussion the principle of exchange of 
stabilities is assumed to hold So much so putting 0 = 0, tiie 
set of Ihuations (5 11 ) and (5 12) becomes 

2v2 ^f'^f 

(D - a ) = — 5- a e 
f 




(D^- a^) 6 = - 


(I.^- a^)w-= - " “ 


B 1C Aa e' 


m 


2 ° i 2 id 


„ yB d*^ 

(])^- a^)0'= w' 


(5 14) 


(5 15) 
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Ihe conespoMing bovindaiy coEditions are 


At 

z 

= 0 

e* 

= w‘ = 0 


At 

z 

. h 
= 1 

0 

=s w =0=Dw 


At 

z 

II 

0 

- 0 ', w w', 1 ^ DO' = D0 

(5 16) 


/ 2 2 1 
(D - 2a )Dw + - 7 - Dw' =0 
eA 

2 “0 

D w s= (Dw - Dw') 

/eAT 

In order to apply the above method, the following four 
functions are defined 


w = (z - 1 (a^+ a^z + z^) 

9 = (z - 1 -|) (b^+ z) (5 17) 

W*= (Cq + z) z 
6'= z(d^+ z) 

lEhese are the ’test' functions and satis^ part of the boundary 

conditions !Ihe constants a . a., b^, and d^ are determined 

0 1' O' o 0 

by the remaining conditions Hieir values are 


^ = (u/2- 

% = (*'2- 


= a, + dA 

h 




'•m 


^0 = 


m i 


(5 18) 
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where 


^ /il /il ^ 


TH- i / 2 h"' h^ / i; 

Ifg - 12 a - ♦=> -^ + -j (-ta-^. 
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= ("V + 't) 

eA d-' 
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ii 


•^'' A /' H h\2n / 2 a 

J \ (z - 1 ” j) (a^ + a^z + z ) 


2 An 


0 = I Ba (V 1 -|)“ 

n=1 

w’= y 0 (c + z)^ z^ 
n=1 

0*= I \ z^(d + z)^ 

it=1 


(5 19) 


As is done in the G-alerkin’s method, these expression are substituted 
in the set of Equations (5 14) and (5 15) Since these are not the 
actual solutions, these will not satis:^ the differential equations 
(5 14) and (5 15) So much so the substitution of these ex:pressions 
into the differential equations will result in residuals E^, r(w) - 
Q(6) = Ep where P,Q are polynomials in z coming from tenas 
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depending on w ard e respectively IJie residuals corresponding 
to each of the equations of the set (5 14) and (5 15) are normalized 
!Ihis will lead to a set of four homogeneous linear equations in the 
coefficients C^, In order that a non-tnvial solution 

exists for these coefficients, an infinite ordered determinant has 
to vanish sucessively taking the first, second etc terms xn 

the expansion, one can determine the critical Rayleigh numbers Infact 
the Rayleigh numbers will satisfy functional relations of the form 


R^ = f ( 3-q t j , c^j a) 

(5.20) 

Rjjj^ = S (S'q# °o» *^0**^0*^* d" * 

From these the critical Rayleigh numbers can be determined 

Ihe evaluations were done taking the first and the first 
two terms of the expansions ®ie calculations are indicated in 
the following Bee use of the procedure outlined in the previous 
paragraph will result in the following system of equations for the 
case vdien the first terms are used 

3jC 


(24S(1) - 2a^ v;, + ^ ° 


(5 21) 




^ 7 *^ + (2u(l) - = 




Y? d' 
m 


(5 22) 
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Y0 d o^d 
m 




m 


5 'd' ‘ 4 '■d 

.2 


d 

i 

*3 




where 


V*, -QiiS( 1 )+«i 23 ( 2 ) + Qi 3 S( 3 ) 


+ I,2S(2) + I„S(3) + 1 , 4 ( 4 ) S(4) 4 S(5) 


•12 

=H„ (1) 4 E, 28(2) 4 S(3) 

Y* - B„tj(l) 4 E ,D(2) 4 0(3) 

4-1 M ‘ 

. ]i + li!.p-.+ ^ 2 P 
E..=V®' 1 d + ^2 M 2 1 d 


11 


•11 


P - 1 - 2(a^+-|^ ) + P 


•15 

2h 

P . = ■+■ "7” 

■^14 1 ^ 


11 


Q1I “ ^ ^13^ ^12 


> <S« - «*14 “13= ' = 


1 2^3. 4x i-1^ ^^12 ,9 

+ P.,2 + ^ ' 


D(l) = / (:Pii-^^ 12 "-^-"i 3 '^ 


1 2\ 1 1 1 - 1.2 ,9 

u(x)==/ (Rii-^ ^12=^ ^ ^ ^ 

0 

„ b o,a ar as glvanittBluatlon(5 18) 
where a^i a^» ®q» %* 0 

(5 21) and (5 22) one btain 


Prom UQuations 



57 


Ej = - (245(1) - aS*^) ( 20 ( 1 ) - 

2 


(5 23) 


R 

m 


= [(Oo.f|)-a^^(|)^lo^(|)^,l(|)3,J 

^2 

^ J ^ ^“1 ■*■ i F divided by 


AaW[! 2!o 2 

'd^ *■ 3* 4 d~52'^ 

d 

VJhen. the first two terms in the expansions are used, the following 
equation is got for the Rayleigh njmber R™ 


R. ^ (-V32 + *^(V32f - 4 V31 V33) / 2 733 


(5 24) 


Here 

731 = (V25 716 - 715 726) (711 722 - 712 721 ) 

732 = -[(721 728 - 727 722) (713716 - 715 714) 

+ (721 718 - 717 722) (725 714 - 713 726) 

+ (727 712 - 711 728) (723 7l6 - 715 724) 

2 

+ (711 718 - 717 712) (723 726 - 724 725)1 a 

733 = (723 714 - 713 724) (717 728 - 727 718) ^ 

VII = 24 S(l) - 2a^ {Q11 S(l) + Q12 s(2) + Q13 S(3) > 

+ {P^^S(l) + I^2®(2) +1^33(3) + P^4S(4) + S(5) } 

712 = 24 {S13 S(l) + 5 S12 3(2) + 15 311 s(3) - 70 PI4 S(4) 
+ 70 3(5 ) } - 4a^ { 315 S(l) +3 Sl4 S(2) +6 SI3 8(3) 
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+ 10 S12 S(4) + 15 311 S(5) + 42 P14 S(6) + 28 S(7) } 

+ {P17 S(l) + S16 S(2) + 315 8(5) + S14 8(4) + 813 8(5) 

+ 812 S(6) + 811 8(7) + 818 S(8) + 8(9) > 

V13 = R11 S(l) + R12 8(2) + 3(3) 

?14 = RI3 S(l ) + R14 8(2) + R15 8(3) + R16 8(4) + S(5) 

V15 = 2U(l) - {R11 U(l) + R12 U(2) + u(3) > 

V16 = R17 U(l ) + R18 U(2) + R19 U(3) - a^{E13 U(l) + E14 u(2) 

+ R15 U(3 ) + R16 U(4 ) + U(5 ) } (5 25 ) 

V17 = P11 U(l) + P12 U(2) + P13 U(3) + P14 U(4) + u(5) 

V18 =; P17 ir(l) + 816 U(2) + Sl5 U{3) + 814 Il(4) 

+ 813 U(5) + 812 U(6) + 811 U(7) + Sl8 U(8) + U(9) 

V21, V22, V23, V24 are got by replacing s(i) by t(i) in VII, V12, 

VI3, V14 and V25, V26 , 727 , 726 are got by replacing U(i) by e(i) 

in Uiuation (5 25) The expression for T(i) and e(i) are 

1 2545 

= / [P17 + SI6 z + 815 z + 814 z + 813 z + Z12 z 

0 

6 7 8,1-1, , , q 

+ 811 z + 818 z + z J z dz , 1 = 1, 

R(i) = / (R13 + Rl4z + R15 z^ + E16 z + z ) z dz, i=1, ,9 

o 

where 

S11 = (6 - 8 P15 + PI6 + 2P1l) 

812 = 2(-2 + 6 P15 - 2P16 + PI8 - 4P1l) 
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S15 = 1 - 8 P15 + 6 Il6 ~ 8 P18 + 12 P11 + P17 

514 = 2(P15 - 2 P16 + 6 P18 - 4 P11 - 2 PI?) 

515 = P16 - 8 P18 + 2 P11 + 6 PI?) 

516 = 2(P18 ~ 2 PI?) 

SI? = 1 

S18 = 2 P14, 


P15 = a^+ 2h/^ 

P16 = (a,+ 2h/d)^ 

/ h h \2 

PI? = (a^+ + -j) 

a 

p,e = {a,, f ) (a„. a J * J) 


r 

R13 = (b^+ h/^y 


h^2 


E14 =» 2(b^+ b/d) - 2(b^+ -j) 


R15 « (b^+ jf - 4(b^+ |) + 1 


R16 = 2(b^+ -f ) •- 2 


R1? = 2 R15 
R18 == 6 Rl6 
R19 =* 12 

I 

Ihe remaining yanables have already been defined 


The critical Raleigh numbers R^ and R^ corresponding to the 
fluid and the fluid m the porous medium have been calculated numer^ 
ical3y for different values of the parameters A, ''f 
Equations (5 23)and (5 24) ’Ihw have been graphically shown in 
Pig 5,2, 5 3 and 5 4 A representative sample of these values are 
Tables (5 1 ), (5 2) , and (5 3) 


also given in 
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Dxsoussxon 

She Gz^idMBl B^leigh numbers were deteimxned for different 
values of the parameters !Ehe values of the paramete2?s were chosen 
to be 0 1, 0 12, 0 13, 0 14, 0 15, 0,16 for 2000, 1000, 600, 

5 00, 400 , 300, 250 , 200, 135, 100 for 1 //a and 0 1, 0 5, 8, 1 0, 

1 5, 2 0 for 3he calculations were also perfomed for 

different values of ranging from 0 1 to 1 1 0 The critical 
values of the Rayleigh number are shown in Tigs 5 2, 5 3 and 
that of 5 4 for some of the above values of the parameters 

Erom the figures one observes the following Ttor R^, the 
critical Raleigh numbers continuously increase with increasing 
values of I/v'A justifying the explanation of Sparrow et al Tbr 

R * the critical values increase quite rapidly with i//A while they 

m^ 

seem to tend to a constant value for the case of R^ Infact one can 
easily determine from Equation (5 23) that increases to infinity 
as A increases to zero (l/A “ ) whenever the parameter 
is a non-negative finite constant Also the critical I^’s are far 
greater than those of R^s which m^jlies that the instability sets 
in the fluid layer faster than in 14ie medium Also the curves for 
larger a^areabove those for smaller ct^ 

Qualitatively, the results are as follows for the different 
values of the parameters The values of the critical nuabers are 
highly dependent on the ratio b/d Tbr example, and R^^^ were 
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fouM to be of thi. order of lo"'* and 10^ for h/d = ii while for 
h/d = 0 1 these were of the order 10^ and 10^ respectively Star 
increasing values of h/d, the critical values decrease in a propor- 
tionate fashion* Wien h/d is small, the saturated porous block is 
quite thin as compared to the fluid Iqyer over it Ohe permeable 
boundary is probably affecting the stability of the fluid l^er only- 
through the slip boundary condition Bae patterns within the porous 
material donot probably liave much influence on the fluid layer above 
Brat explains the closeness of to critical Rayleigh numbers in 
the B'^nard problem Hrwever, if h/d is large, the flow behavior 
within the porous material wLll have significant effect on the 
stability of the fluid layer Waich will be very thin compared to the 
block Bor large h/d, a more careful trea-fcment within the porous 
material will be desirable for ary concluBive results As "tiie 
principal aim of the present investigation is to know the effect of 
the slip condition on Benard problem, more attention has been paid 
to -fehe case when h/d is small 

Ihe results of the 'second iteration’ show that the convergence 
of the Galerkin method used is fairly good One sees that for certain 
range of 1 //a , the critical numbers for the 'first' and 'second' 
iteration differ only is the unit place One then can anticipate 
the result to converge in a wider range of 1 //a if more and more terms 
are taken in the expansion of the Oaleakin method# 
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Oiie critical values of 10^ x R for different a ’s are also 

m 0 

close As such only the values for a^= 0 i6 and 0 i have been 
given in the [Cable 5 3 
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Fig. 5.1 BeNARD proslem with porous boundary 





GriAITEH - VI 


SOME COMM0ITS (® THE THESIS T^IG 


The discussion about the problems dealt with In the 
thesis has been mostly in a qualitative manner The reason 
18 simple it is because of the non -availability of the 
exact boundary condition(s) at the interface of a porous 
material for developing flows The proposed condition is 
rather intuitive and no rigorous proof for the ccndition 
has yet been got Possibly some light would be thrown on 
the condition with some experiments of developing flows on 
porous material Once the correct boundary condition(s) 

IS available, a numerical integration of the full Ifavler-Stoifces 
equations would give a clearer picture of the entry flow in 
a porous walled channel 3his would lead to accessory problems 
like the convergence of the method used, stability, and the 
computer memory needed Again one ims to discuss the existence 
and uniqueness of the solution for the PDE 

Agam for the problem of flow over the porous flat 
plate, no proof of existence and uniqueness af solution has 
been found as yet A proof on the lanes referred to by 
ITickel (1973) in his review paper may be possible The second 
solution for the problem has been attempted m this ligjit 
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Ibr the Benard problem, the principle of exchange of stability 
has not been, proved But for the two mediums - fluid and porous - 
the proof of such a phenomenon is available of course using the no-slip 
conditions 
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